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Simulations of momentum feedback by black hole winds 



o 
o 

<N 
> 

o 



o 

u 

Ok 

6 

h 



> 
m 

<N 



o 



Sergei Nayakshin 1 ' 2 and Chris Power 1 

1 Department of Physics & Astronomy, University of Leicester, Leicester, LEI 7RH, UK 

2 Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Strasse 1, D-85741, Garching, Germany 



Accepted 2008 ?? ??. Received 2008 ?? ??; in original form 2008 05 ' 



ABSTRACT 

The observed super-massive black hole (SMBH) mass - galaxy velocity dispersion 
(Mbh — cr) correlation may be established when winds/outflows from the SMBH drive 
gas out of the potential wells of classical bulges. Here we present numerical simulations 
of this process in a static isothermal potential. Simple spherically symmetric models of 
SMBH feedback at the Eddington luminosity can successfully explain the Mbh — cr and 
nuclear cluster mass Mnc — cr correlations, as well as why larger bulges host SMBHs 
while smaller ones host nuclear star clusters. However these models do not specify how 
SMBHs feed on infalling gas whilst simultaneously producing feedback that drives gas 
out of the galaxy. 

More complex models with rotation and/or anisotropic feedback allow SMBHs to 
feed via a disc or regions not exposed to SMBH winds, but in these more realistic 
cases it is not clear why a robust Mbh — cr relation should be established. In fact, some 
of the model predictions contradict observations. For example, an isotropic SMBH 
wind impacting on a disc (rather than a shell) of aspect ratio H/R <C 1 requires the 
SMBH mass to be larger by a factor ~ R/H, which is opposite to what is observed. 
We conclude that understanding how a SMBH feeds is as important a piece of the 
puzzle as understanding how its feedback affects its host galaxy. 

Finally, we note that in aspherical cases the SMBH outflows induce differential 
motions in the bulge. This may pump turbulence that is known to hinder star forma- 
tion in star forming regions. SMBH feedback thus may not only drive gas out of the 
bulge but also reduce the fraction of gas turned into stars. 
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1 INTRODUCTION 

It is believed that the centres of most galaxies con- 
tain super-massive black holes (SMBHs) whose mass 
Mbh correlate s with the velocity disp e rsion a of the 
host galaxy (|Ferrarese & Merrittl l200d : iGebhardt et al.l 
l200d : iTremaine et all l2002h ~ Similarly, there is a cor- 
relation between Mbh and the mass of the bulge, 
Mbuiee for la r ge SMBH masses (|Magorrian et al.1 Il998l : 
Har ing fe Rixl 120041 : iGiiltekin et al. I l20Q9h . Observations 
also suggest that the masses of nuclear star clusters 
(NC) (1O 5 M < M NC £ 10 8 M o ) cor relate with the prop- 
erties of their host dw arf ellipticals ([Ferrarese et al 1 l200d : 
IWehner fe Ha rris 2006) in a manner that is analogous to 
the one between SMBHs and their host ellipticals. 

These empirical relations can be explained in a very 
natural way if the growth of host galaxies and their central 
SMB Hs or NCs are linked by feedback. This was first pointed 
out bv lSilk Reesl (|l998h , who highlighted the potential im- 
portance of SMBH heating and outflows before any robust 
observational evidence for such a link had been found. Sub- 



sequently, iFabianl (|l999h argued that radiation pressure act- 
ing on cold gaseous cl ouds in the bulge could give rise to the 
observed correlations. iKing fe Poun ds (2003) argued for the 
existence sub-relativistic outflo ws from the very c entral re- 
gions of AGN, which prompted iKingl fcOOl l2005h to study 
the motion of a shell of gas swept up by a wind/outflow 
from a centra l black hole in a galactic isothermal dark mat- 
ter potential. IKingl (|2QQ5h demonstrated that the shell will 
be expelled from the potential provided the black hole mass 
exceeds a critical value that, as a function of cr, turns out to 
be close to the observed Mbh -a relation. 

The main result of kind (|2Q0l l2QQ5h can be deduced 
using a simpler order of magnitude "weight argument" . Ac- 
cording to this argument, the SMBH luminosity is assumed 
to be limited by the Eddington valu e. Radiation pressure 
drives a wind. IKing Pounds] ([2003) argue that wind ve- 
locity is comparable to the escape velocity from the inner 
accretion disc, e.g., v ~ 0.1c. The momentum outflow rate 
is assumed to be 
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here k is the electron scattering opacity and Mbh is the 
SMBH mass. This result is natural to the order of mag- 
nitude: LEdd/c is the radiation momentum flux which is 
presumably passed to the wind as radiation accelerates the 
outflow (cf. iKing fe Pounds! l2003h . The argument also as- 
sumes that the black hole wind is optically thin at the 
point of interaction with the ambient gas. Because the cool- 
ing time of the shocke d gas is short on scales appropri- 
ate for observed bulges (King 2003L l2005h , the bulk energy 
of the outflow is thermalised and quickly radiated away. 
It is then only the momentum push (equation [T]) of the 
outflow on the ambient gas that is important since it is 
this that produce s the outward force on the gas (as in the 
earlier model by Fabian 1999). The weight of the gas is 
W(R) = GM(R)[M tota i(R)]/R 2 , where M gas (R) is the en- 
closed gas mass at radius R and M to tai(-R) is the total en- 
closed mass including dark matter. For an isothermal po- 
tential, Mg as (R) and Mtotai(-R) are proportional to R, so 
the result is 



W : 



G 



(2) 



Here f 9 is the baryonic fraction and a 2 = GM to tai(R)/2R 
is the velocity dispersion in the bulge. By requiring that 
momentum output produced by the black hole just balances 
the weight of the gas, it follows that 
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(3) 



which is consistent with the observed the Mbh-0" relation. 

To order of magnitude, relation [2] should hold fo r any 
poten ti al at the virial radius. Therefore the model by IKing 
(12001 l2005h appears to be a promising explanation of 
the observed correlations between SMBHs and their host 
g alaxies . How e ver , a dditional complications, not considered 
m iKind (I2QQ3I . [2QQ5h . may be important. Amongst these are 
(1) the self-gravity of the gas, because gas that accumulates 
at the centre of the potential may begin to dominate it and 
therefore alter the result; (2) finite angular momentum of 
the gas may lead to formation of a disc; (3) collimated and 
variable outflows. 

The goal of this paper is to v erify t hat th e predictions 
of the analytical model of IKind (I200I l20Q5h hold and to 
explore more realistic settings, of the kind just described. 
Note that we concentrate on the early stages of galaxy evo- 
lution, when we would expect galaxies to be gas-rich and 
the baryonic mass within dark matter haloes is dominated 
by gas rather than stars. At later times, when the SMBH is 
likely to be less luminous and the galaxy is less gas-rich, we 
would expect different forms of feedback to become impor- 
tant, such as relativistic jets an d the associated radio bub- 
bles (e.g., IChurazov et al. I l2002h or pre-heating due to the 
inverse Compton effect (e.g. JSazonov et alj|2005h . However, 
we do not include these forms of feedback in our current 
simulations. 



2 FIXED BH MASS RUNS 
2.1 Numerical method 

iNavakshin et all (|2009aT ) developed a new method for ra- 



diation transfer in SPH based on Monte Carlo packets, 
which can also be used to simulate winds in the momentum- 
conserving phase. As discussed in the Introduction, one can 
consider only the momentum transfer from the outflow to 
the ambient gas in this regime. In addition, for the prob- 
lem at hand, the typical velocity of the ambient gas is of 
the order the velocity dispersion in the bulge, which is less 
than a few 100 km/sec. In contrast, the black hole outflow 
is much faster, v ~ 0.1c, which implies that the mass of the 
black hole wind is negligibly small compared with that of 
the gas in the bulge at the point there the latter is driven 
away. Hence we can neglect the mass of the black hole wind 
in comparison to that of the ambient medium. 

Accordingly, we use massless particles moving with ve- 
locity v = 0.1c to simulate the black hole wind. The wind 
particles propagate in straight lines until they encounter 
one or more SPH particles at which point they transfer 
their momentum to the particles. For mor e information on 
the me thod and validation of the code see INavakshin et all 
(l2009aT ). especially their §3. 

For spherically symmetric tests, we use 35,000 SPH 
particles. The wind particles carry momentum p w ind = 
0.1ra sp hcr each. This satisfies the requirement p w ind <C p sp h, 
where p sp h is the typical SPH particle momentum (~ m sp h<J 
here), and ensures that Poisson noise from our Monte Carlo 
scheme is small enough not to compromise our results (see 
Nayakshin et al. 2009a). The particles are emitted by the 
black hole at the rate 

jV wind = , (4 ) 

CPwind 

where L^dd is the Eddington luminosity for the black hole. 
They are ejected isotropically unless stated otherwise. 

We perform our simulations in a static singular isother- 
mal sphere potential with velocity dispersion a = 147 km 
s~ . The units of mass and distance are Mu 1O 1O M 
an d Ru — 1 kpc respectively; the unit velocity is Vu = 
GMu I Ru = V%cr. Equation [3] yields the expected value 
for the critical mass at which gas should be expelled; 
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1.1 x 1O 8 M 



(5) 



To check this prediction, we first limit our simulations to 
spherically symmetric initial conditions. We work with a fi- 
nite extent shell since that allows us to explore interesting 
regions of parameter space. The initial radial distribution 
of gas follows the distribution of dark matter, assumed to 
dominate the potential: p(R) oc R~ 2 for R in < R < R Q 



For a ll the tests in this paper, Ri n = R out /2. As in IKind 
(2003|, 120051 ). the normalisation of gas density, and hence 
the mass of the shell, is set by the cosmological baryon mass 
fraction, f g0 = p g /ptotai — 0.16, where ptotai = Pdm + p g 
is the total mass density, and p g and />dm are the gas and 
dark matter densities, respectively. For the isothermal po- 
tential, ptotal (R) = (j 2 /2tvGR 2 . Because our focus is momen- 
tum feedback, we adopt the isothermal equation of state for 
the gas, assuming T = 3 x 10 5 K. This is lower than the 
virial temperature, T v i r ~ 10 6 K, for the potential we are 
using, preventing small scale gas fragmentation. The latter 
process would lead to star formation and feedback from star 
formation, in both heating and momentum-driven forms. We 
plan to include these processes in our future work, for now 
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neglecting them in comparison with the SMBH feedback al- 
together. 

We set Race — 0.5 as the inner boundary of our 
computational domain. SPH particles that cross inside this 
inner boundary are removed. Note however that this scale is 
still too large to associate it with SMBH accretion directly, 
and there is no model-independent parameter-free way of 
connecting this "sunken" mass with SMBH accretion and 
feedback. For simplicity we take the pragmatic approach 
in this paper and consider only tests in which the black 
hole feedback is either constant or increases at a rate fixed 
by the Eddington limit. This allows us to concentrate on 
gas dynamics without worrying about the complexities of 
the non-linear and model-dependent accretion-feedback 
link. The "accreted" gaseous mass is thus considered to be 
mainly used in star formation inside R acc and the feedback 
from that is neglected. 

Before we begin the discussion of results, we show in 
Figure [T] two typical snapshots from the spherically sym- 
metric calculations. These show gas column densities in an 
angle-slice projection. Specifically, the gas column density 
shown in the Figure [T] and in all similar figures below, is 
calculated by 

/z(x,y) 
p(x,y)dz , (6) 
-z(x,y) 

where the limits of the integration are given by z(x,y) = 
rtan£, and r = (x 2 + y 2 ) 1 ^ 2 . The angle £ is chosen to be 
tan£ = 1/4 for this paper. This projection method is con- 
venient as it permits an unobstructed view into the central 
regions where the black hole resides. Had we used a constant 
thickness z(x,y) = const projection method, then either the 
innermost region would not have been resolved sufficiently 
or an insufficient number of the SPH particles in the out- 
ermost regions would have been sampled for a statistically 
meaningful figure. 

The tests shown in Figure [T] are described in g23] (left 
panel) and £12.41 fright panel). Briefly, the former has Mbh = 
2 x 10 8 M , whereas the latter has M bh = 4 x 10 8 M . In 
both cases the gas is falling in with radial velocity v r = — 1 
(in code units) initially. In the lighter SMBH case, the black 
hole momentum outflow is insufficient to reverse the shell's 
infall and it engulfs the black hole. In the right panel, on 
the other hand, the initially thick shell is first compressed 
to a thinner shell by the opposing actions of the black hole 
outflow and the inward inertia of the outer layers, and then 
expelled to infinity. 



2.2 Initially static shell tests 

We begin with tests in which gas is initially at rest (even 
though this is unrealistic for the chosen temperature). The 
initial outer radius of the shell is R on t = 40 kpc for these 
runs. Figure [2] shows results of three such tests, where 
we vary the black hole mass from Mbh = 5 x 10 7 M to 
Mbh = 2 x 10 8 M by factors of 2. We define a mean ra- 
dius and velocity of the gas by averaging over all the SPH 
particles in a simulation. Defined in this way, the "shell" 
radius is plotted in units of the initial mean radius, whereas 
the velocity is shown in code units. It is apparent from Fig- 
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Figure 2. Mean radius (upper panel) and radial velocity (lower 
panel) of the gas as a function of time for three different black 
hole masses: M bh = 2 x 1O 8 M , solid,M bh = 1O 8 M , dotted, 
and Mbh = 5 x 10 7 Mq, dashed. The radius is plotted in units of 
the initial mean radius, whereas velocity is shown in code units 
(V2a). See 322] for more detail. 

ure[2]that the lower mass black holes (Mbh ^ 1O 8 M ) are 
unable to expel the shell of gas, whereas the more massive 
black hole (M bh < 2 x 10 8 M ) is able to expel the shell. We 
therefore estimate that the critical black hole mass for this 
potential is around 1.5 x 10 8 M , in a very good agreement 
with the estimate in equation [5] 

2.3 Initially infalling shell, M bh = 2 x 10 8 M 

This simulation is identical to that presented in £12.21 with 
Mbh = 2 x 10 8 M , except the radial velocity of gas is now 
non-zero, with v r = —1. The additional initial inertia of the 
shell suggests that the critical M a mass should be higher in 
this case, which turns out to be true. 

Figure [3] shows the mean shell radius, i? s h, in units of 
its initial value (dash-dotted green curve) as a function of 
time. The inward motion of the shell cannot be prevented in 
this simulation, as is clearly indicated by the monotonically 
decreasing mean radius. The solid curve shows the mean ra- 
dial velocity of the gas scaled to its initial value. The shell is 
slightly decelerated initially by the momentum of the SMBH 
outflow, i.e., |i7 r | decreases until t ~ 40 Myrs. However, at 
t ~ 40 Myrs the shell starts to accelerate and falls "onto" 
the SMBH, completely enveloping it. To follow the evolu- 
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Figure 1. Typical results of the spherically symmetric tests. Left panel: shell collapsing on the Mbh = 2 x 10 8 Mq black hole ( 32.3) 1. 
Right panel: same shell but the black hole mass is M b h = 4 x 10 8 Mq, sufficient to drive the shell out Both snapshots are done 

at time t = 70 Myrs. See 32. II for description of the figure. 



tion of the shell after this point would require much greater 
numerical resolution and, more importantly, the physics of 
star formation. 

The additional inward acceleration of the shell at t ~ 40 
can be understood if one considers the self-gravity of the gas. 
We calculate the velocity dispersion corrected for the shell's 
self- gravity: 



■I 1 



CD 

= ^ [Mum(R) + M Bh ] = (1 - f g0 ) ol + ^ , (7) | 

where R = R s h- The second term increases for a constant ^ 
M s h and a decreasing R s h. It is convenient to define a "run- 
ning" value of the expected M a value based on the mean 
values of cr s h and R s h using equation [3] Further, we use the 
dimensionless variable ma' 



M(r(R s h, (7sh) _ fg ( CT s h 



M BH =2x10 8 M sun , R out =40kpc 
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where M a o — M a (Ro, ao) is the initial value of the M a mass 
as a function of the initial values of shell radius Ro and veloc- 
ity dispersion ao , and f g and f g is the current mass fraction 
of gas inside radius R. The function is plotted in Figure [3] 
with the dashed curve for the present simulation. It becomes 
immediately apparent that the SMBH will find it harder to 
retard the shell once the shell has fallen in sufficiently deep 
into the potential and the gas starts to dominate the gravi- 
tational potential. If SMBH mass is below the critical mass 
then a runaway radial contraction of the shell occurs in this 
simple model. 



Figure 3. Mean radial velocity (solid) and radius (dash-dotted) 
of the gaseous shell as a function of time for Mbh = 2 x 10 8 Mq. 
The dashed curve shows the expected dimensionless mass, as 
defined by the equation [8] The simulation is described in £12.31 

2.4 Initially infalling shell, M bh = 4 x 10 8 M 

Figure [4] shows the result of exactly the same calculation 
but now with a SMBH mass twice as large, i.e., Mbh = 
4 x 10 8 Mq. The behaviour of the shell in this case is radi- 
cally different. The average radius of the shell decreases only 
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Figure 4. As in Figure [3]but now the black hole mass is a factor 
of two larger ( ^2.4|) . The gaseous shell is now completely expelled 
from the potential well to infinity. 

during the first 40 Myrs, during which time the outermost 
layers of the shell continue to fall in. However, at t ~ 40 
Myrs the shell briefly decelerates before quickly accelerat- 
ing nearly linearly (v r oc t) out of the potential well. 

It is worth noting how distinct the gas dynamics be- 
tween the two cases illustrated in Figures [3] and U re- 
spectively. This may be somewhat surprising given that the 
SMBH mass and hence the outward acceleration differs be- 
tween the cases only by a factor of two. However, the key 
difference is that in the first case the shell fell in far enough 
to become self- gravitating, with the self-gravity corrected 
a increasing rapidly, increasing further the inward p ull of 
gravity. This effect, missing in the analytical theory of iKind 
(2003j, I2005I ). might actually increase the M a value above 
that given by equation [3] 

Comparing the two initially infalling shell tests with 
those for an initially static shell ( fl2.2[L we note that the 
critical black hole mass is increased by a factor of two or 
so to M a « 3 x 10 8 M©. This is natural as the shell needs 
to be accelerated to a positive velocity v r ^ 1 to be driven 
out of the potential, and that is comparable to the initial 
negative velocity of the infalling tests. The required outward 
acceleration should then be roughly doubled compared with 
the initially static tests. This shows that variations in the 
initial radial velocity of the shell are unlikely to change the 
critical black hole mass by more than a factor of a few. 

2.5 Scale- free nature of fixed BH mass solutions 

Equation [3] shows that the critical M a mass is independent 
of the initial location of the shell, i?o, as long as the ini- 
tial gas fraction f g o = M gas (i?o) /Mtot(Ro) does not vary 
with radius. We now show that the same is true for time- 
dependent solutions. This can be seen from the equation of 
motion for gas in the one zone approximation, which is quite 
reasonable for a thin shell. We have 



Figure 5. As in Figure 0] but now the initial shell radius is 10 
kpc rather than 40 kpc. Note that the result is identical to Figure 
[4] save for the time axis shrinking by factor of 4. This confirms 
the scaling properties of the model discussed in £12.51 



M sh ^ = ^ - |L [Mum(R) + M sh ] M sh , (9) 
at C K 

where v r — dR/dt. Let us assume R(t) = Rog(^), where 
£ = t/to, to = Ro/a. Obviously, we require g(0) = 1. We 
can now re- write the above equation in this form: 

R» #g = ^Edd _ 2a 2 (l-/ g0 ) _ GM sh 

tl de cM sh R g R 2 Q g 2 ' 1 } 

Only the dimensionless time variable £ remains, implying a 
scale- free solution. 

This conclusion is confirmed in Figure [5] which shows 
the result of a simulation identical to the one presented in 
Figure 2] except that the initial outer shell radius is .Rout = 
10 kpc rather than Rout = 40 kpc. The shell mass is also 
reduced by a factor of four to keep f g o the same. It is readily 
apparent that the two figures are identical, with the only 
difference evident in the time coo rdinat e. 

This scale-free nature of the iKind (|2003l . l20Q5h model 
in the isothermal potential has interesting observational im- 
plications. The model predicts a slope and normalisation for 
the Mbh — cr relation that is consistent with the observed 
ones, assuming only one free parameter / g o fixed at the cos- 
mological baryon fraction; however the overall scale of the 
syste m needs to be set by additional physics. IKind {2003, 
I2005T ) points out that this can be the requirement that the 
black hole outflow's cooling time is short compared with the 
local dynamical time, which fixes the ratio of the black hole 
to bulge mass to ~ 10 ~ 3 , in a good agree ment with the 
observational data fe.g. jHaring &; Rixll2004h . 

In principle, other factors can be affecting the black hole 
to bulge mass relation, such as the finite extent of astrophys- 
ical potential wells, or the amount of angular momentum in 
the gas. We shall now consider the first issue in relation to 
the finite maximum growth rate of the black holes. 
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EDDINGTON RATE 
MODELS 



LIMITED GROWTH 



In 32] the black hole mass was assumed fixed for simplicity. 
In this section we allow the SMBH to grow at its Eddington 
accretion rate. In this case the SMBH mass increases by the 
factor of e in the Salpeter time, 



tSalp 



M bh 
M Ed d 



: 4.5 x 10 7 yrs 



(11) 



where MEdd = L^dd/ec is the Eddington accretion rate, 
and the radiative efficiency e is set to 0.1. 

The initial value of the outer radius, R ou t, defines an- 
other important time scale of the problem - the dynamical 
time, 



^dyn 



Rout 

~Vu 



(12) 



where Vu = y GM/R w 208 km/sec is the velocity unit for 
the potential. 

As emphasised by iNavakshin et al.l (|2QQ9bh , there are 
two distinctly different regimes. If tsaip <C tdyn, then the 
SMBH can grow arbitrarily quickly if provided with enough 
fuel. If there is a SMBH feeding- feed back l i nk th at can limit 
the SMBH mass, such as in the lKind (|200l l2005h model, the 
latter then grows to the appropriate M a mass and remains 
there. If tsaip ^> tdyn, then the SMBH is unable to grow 
sufficiently quickly to reach its maximum (i.e. M a \ even if 
it is provided with ample fuel during the dynamical time of 
the system. We now present two simulations that explore 
these two regimes. 
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Figure 6. Gas mean velocity (solid), scaled mean radius (dot- 
dashed), scaled M a (dashed) and the ratio of the SMBH mass 
to the initial M a = M = 3 x 10 8 M mass (dotted) for an 
Eddington limited simulation ( 33.1)1 . The initial SMBH mass is 
10 8 M©, a third of the M a mass. However, due to rapid accretion 
the SMBH catches up with the growing (contracting) bulge and 
expels the gas when M b h « 1.3 x 10 9 M©. 



3.1 SMBH growth in a "large" bulge 

The initial condition used in this test is same as in fc|2.31 ex- 
cept that the initial SMBH mass is smaller, Mbh = 10 8 M©. 
As we found in fc|2.4l for this initial condition, the M a mass 
is about M = 3 x 10 8 M©. The black hole thus needs to 
increase its mass by about a factor of 3 to drive the shell 
out. The dynamical time of the shell is tdyn ~ 200 Myrs, 
which gives the SMBH plenty of time to grow. 

Figure [6] shows the mean radial velocity of gas (solid 
curve), the mean radius of the shell (dash-dotted), the self- 
gravity corrected velocity dispersion of the gas (dashed) and 
finally, the dotted curve shows the ratio of the black hole 
mass to the initial M a = Mo mass. 

The early phase of gas dynamics is quite similar to the 
underweight fixed mass case considered in fc|2.3l The shell is 
contracting and the gas velocity dispersion grows with time, 
increasing the M a value (see the dashed curve in the figure) . 
However, the SMBH grows even faster. At around 90 Myrs 
the gas suddenly gets decelerated and then accelerated to 
positive velocity. 

While this i s as expected based o n simple analyti- 
cal expectations ([Navakshin et al.ll200 9b). the SMBH mass 
at the time when gas velocity becomes positive is about 
1.3 x 10 9 M©, about a factor of 4 higher than the initial 
configuration value, M a — 3 x 10 8 M©. At the same time, 
the shell radius is much smaller at that moment than the 
initial value, increasing the self- gravity corrected M a value 
by a factor of about 2.5 to ~ 7.5 x 10 8 M©. If our mod- 
els included star formation and if a good fraction of gas 



was turned into stars then the resulting bulge, assuming the 
stars remain bound as the gas is blown away, would satisfy 
the Mbh — cr relation within a factor of two. However the 
bulge velocity dispersion would be higher than that of the 
underlying isothermal potential value. 



3.2 SMBH growth in a "small" bulge 

We now repeat the run of fl3.1l but shrinking the shell's outer 
radius by a factor of four to R out = 10 kpc and increasing 
the SMBH's initial mass to M bh = 2 x 10 8 M©. This is about 
two thirds of the initial M a mass, and hence the black hole 
needs to increase in mass by only a small fraction to reverse 
the inflow of gas. However, as Figure [7] demonstrates, the 
black hole's growth is too slow for this configuration of gas. 
As SMBH mass grows, so does the required m a mass, since 
the shell contracts. In fact, when the shell's mass exceeds 
the local dark matter mass the shell becomes self- gravitating 
and the increase in m a accelerates, leaving no chance for the 
SMBH to catch up. 

The results of these experiments confirm that the 
Salpeter time should be sufficiently short compared to the 
the dynamical time of the system in order for the SMBH 
feedback to compensate for the shell contraction. However, 
if star formation had been included, removing some of the 
available gas and producing its own feedback, this require- 
ment could be relaxed somewhat. 
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Figure 7. Similar to the simulation shown in Figure [6] but with 
a smaller initial radius of R ou t = 10 kpc, and a larger initial 
SMBH mass of 2 x 1O 8 M . Note that the feedback from the 
growing SMBH is unable to prevent the contraction of the shell 
in this simulation. 



4 FEEDBACK ON GAS WITH ANGULAR 
MOMENTUM 

So far we have neglected the angular momentum of the gas, 
instead looking at problems in which the motion of the shell 
is purely radial. However, we expect the gas to have non-zero 
angular momentum and this will be as important for the im- 
pact of the SMBH feedback, if not more so, than how the 
the gas is distributed or what its thermodynamical proper- 
ties are. If all the gas in the bulge possesses non-zero angular 
momentum with respect to the black hole, then its feeding 
must proceed through a disc. As efficiency of disc accretion 
is not well understood in the presence of star formation in 
mass ive cold discs (e.g., iGoodmanl 120031 ; iNavakshin et all 
l2007h . the exact distribution of gas in the the angular mo- 
mentum space prior to disc formation is very important. 
Furthermore, on larger scales the angular momentum of gas 
determines the resulting scale of the galactic disc, and per- 
haps the bulge. 

The parameter space of this problem is very large, and 
so we can only begin to scratch the surface of the prob- 
lem here. We choose to study a case of the same initially 
spherical shell infalling with v r = — 1 at t = 0, R out = 40 
kpc (e.g., 32. 3[) , but now rotating around the z-axis. In the 
spherical coordinates, in which z — rcos#, x — rsin#cos0, 
y = r sin 6 sin 0, the rotation law is given by 



Vcji = Vrot sin , 



(13) 



where is the ^-component of velocity, and v ro t — 0.3. 
Thus material near the pole rotates slowly while material 
near the equator (z — 0) rotates at the maximum for the 
test. We expect these tests to be geometrically more com- 
plicated than the spherical shell tests. Therefore, we signif- 
icantly improve the resolution of these simulations by em- 



4.1 Gas dynamics without feedback 

We begin by studying the dynamics of the gas in the absence 
of feedback from the black hole. Figure [8] shows two snap- 
shots of the simulation viewed edge-on. We also show veloc- 
ity vectors of the gas. In this figure and in all others showing 
velocity information, the normalisation of velocity vectors is 
not to unit velocity but to the largest velocity found in the 
snapshot. This is convenient as a subset of gas can have a 
velocity much larger than unity in the more complex cases 
considered below, and its velocity vectors then clutter the 
figures. 

As expected, a gaseous disc forms with an outer radius 
of about R ~ Rout/ 3 ~ 13, appropriate for the rotation ve- 
locity equal to about a third of the circular velocity in the 
isothermal potential. The inner radius of the disc goes all 
the way to our inner boundary of the simulation domain, 
R = 0.5 kpc, since the regions close to the polar axis have 
small angular momentum. By late times about 4% of the 
gas disappears inside the inner boundary. The disc under- 
goes radial oscillations for many orbits before the gas settles 
on circular orbits. Signatures of these oscillations can be ob- 
served in the right panel of Figure [8] where the disc vertical 
structure is not completely relaxed. In particular, one can 
see some low density regions having a small positive radial 
velocity in the corners of the Figure. 

Due to the high gas temperature selected for the simula- 
tions, T = 3 x 10 5 K, the vertical pressure scale-height of the 
final disc is relatively large, H/R ~ 0.3. Realistic discs are 
expected to be cooler, and thus could be geometrically thin- 
ner and denser. However star formation in these discs and es- 
pecially feedback due to star formation may make such discs 
vertically more extended than one would get from simple es- 
timates based on their temperature only. In any case, our 
deliberate choice of a rather geometrically thick disc implies 
that our main conclusions on the role of angular momen- 
tum, to be drawn below, can only be strengthened if discs 
are thinner and denser. 



4.2 Feedback and disc formation 

We now turn to simulations with black hole winds. In or- 
der to isolate the main effects in the interplay between the 
feedback and angular momentum, we consider the simplest 
case of a constant and isotropic feedback. The momentum 
flux from the black hole is fixed at the value appropriate for 
a black hole mass of M bh = 2 x 10 8 M , as in 3231 

Figure [9] shows two snapshots from the simulation. The 
left panel of Figure [9] demonstrates a rather obvious point 
- that inflow is slower in the presence of feedback than in 
its absence (Figure [5}; in fact, the inner region is devoid 
of gas. As in the spherically symmetric feedback cases, a 
shell is formed. The shell is not spherical. The equatorial 
regions are held back by the centrifugal force more than 
the polar regions. Whilst this is occurring, gas also starts to 
accumulate in the midplane, just like it does when there is no 
feedback. Because there is a competition between feedback 
and centrifugal forces on the one hand, and gravity of the 
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Figure 8. Column density in the edge-on projection of the collapsing shell for simulation without feedback. Left: time t = 95 Myrs. 
Right: t = 230 Myrs. The final relaxed disc is hydrostatically supported at H/R ~ 0.3. Note the high density inner disc. About 4% of 
the shell has been accreted by the inner boundary at R = 0.5 kpc. 



other, a ring - rather than a disc - is seen to form just 
beyond the shock front in the right panel of Figure [9] Also 
note that the cavity is not perfectly azimuthally symmetric 
at later times (right panel), due to the development of a 
combination of instabilities and the self-gravity of the gas. 

Figure [10] shows two snapshots at a later time. A bipo- 
lar initially bubble-like structure develops around but not 
exactly on the symmetry axis. The bubble quickly bursts 
and a cavity is opened along the symmetry axis. In that re- 
gion some gas is driven sideways, by the centrifugal force, 
and some to infinity by the black hole wind. At later times 
the central region R ^ 10 — 15 is devoid of gas except for 
the black hole wind. Gas on the interface between the void 
with the higher density gas is continuously stripped away, 
feeding the (secondary) outflow. 

The centrifugal force rarefies the region of the shell near 
the poles, and the black hole winds then evacuate that region 
by pushing the gas away. In the equatorial plane, on the 
other hand, the disc becomes too dense for the black hole 
feedback to have much of an effect. The disc manages to 
shadow its interior from the black hole influence. 

Analysing the net effect of the finite angular momentum 
of the shell on the feedback-shell interaction, we note that 
the polar regions can now be expelled more easily. If these 
regions are to ultimately provide the SMBH with fuel then 
this would reduce the required M a mass. However, it is now 
much more difficult for the SMBH to affect the self-shielding 
disc. The mean density in the disc is ~ R/H times higher 
than the density for the same mass spread in a spherical 
volume. If it is the disc that provides the SMBH with fuel, 
then one would predict a value for the M a mass higher than 



the equation [3] by at least the factor of R/H. Indeed, in case 
of isotropic feedback, the fraction of SMBH wind intercepted 
by the disc is ~ H/R. If feedback is anisotropic and beamed 
away from the disc, as would seem natural for a "grand 
design" disc, i.e., one extending from galaxy scales of kpc to 
sub-pc scales, then there is a further inefficiency in feedback 
delivery to the disc. This appears to contradict the recent 
observations of "pseudo bulges", e.g. bulges that are disc 
or bar-like. The SMBH masses in such bulges appear to be 
smaller rather than larger th an those in classical bulges at 
the same velocity dispersion (|Hull2009h . 

In fact it is hard to see why there would exist any Mbh — 
a relation at all if SMBH were fed by large scale discs that 
are immune to the SMBH feedback. 



5 COLLIMATED FEEDBACK 

There are good reasons to believe that feedback may be col- 
limated rather than isotropic. Unlike stars that radiate due 
to internal energy sources, black holes need a continuous fuel 
supply in the form of accreting gas. Thus there must be re- 
gions where feedback is relatively inefficient and net inflow 
dominates. Most lik ely such inflow takes the fo rm of a stan- 
dard accretion disc (jShakura fe Sunvaevl[l973h on scales of 
mpc and less. The outflow is then probably stronger along 
the axis of symmetry of the disc. In addition, the putative 
molecular torus is probabl y quite a massive and ge omet- 
rically thick structure (e.g., iKrolik fe B egelmar] [l988h that 
may restrict the black hole outflow to the direction perpen- 
dicular to the midplane of the torus (if one can be defined). 
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Figure 9. As in Figure [8]but with feedback, and at times t = 95 Myrs (left panel) and t = 130 Myrs (right panel). 
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Figure 10. As in Figure [9] but for later times. Left: time t = 145 Myrs. Right: t = 275 Myrs. 
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Note that the inner accretion flow does not have to be co- 
aligned w ith the torus, esp ecially if the black hole is rapidly 
spinning ([King et al.ll2005l ) . If direction of the flow fluctuates 
rapidly compared with the time scales of bulge formation, 
~ 10 8 years, then we arrive at possibly quite a complex pic- 
ture of black hole outflows - non- stationary, non-spherical 
and with a fluctuating direction. 

How is the momen tum feedback pi cture developed 
for the spherical model ([Kind 120031 l2005h modified when 
these collimation effects are taken into account? The an- 
swer is probably very complicated and depends on the set- 
ting/environment. Here we present simulations for two rel- 
atively simple situations, designed to motivate further dis- 
cussion rather than give final answers. 

5.1 Outflow misaligned with the disc 

Our first simulation uses a fixed black hole mass M^h = 
2 x 10 8 M© for simplicity, and the same rotating shell initial 
condition used in the simulations presented in 34.21 How- 
ever, rather than assuming isotropic feedback, the black 
hole feedback is now uniformly distributed within the an- 
gle 0.7 < | cos 0' | < 1 around the symmetry axis z' (see 
below). Note that within this conical region the momentum 
flux density carried by the black hole wind is 1/0.3 « 3 times 
higher than in the isotropic case studied in 34.21 

We anticipate that the effect of collimated feedback 
when the axis of symmetry coincides with the z-axis should 
be qualitatively similar to the effect of isotropic feedback 
on gas with non-zero angular momentum, of the kind pre- 
sented in ( 34. 2[) . Therefore we consider a less trivial case in 
which the direction of the axis of symmetry of the outflow, 
z 1 ', is inclined by angle 7r/4 from the z-axis defined by the 
direction of the angular momentum vector of the shell. For 
convenience of presentation, we choose the angular momen- 
tum vector of the shell and the outflow symmetry axis to lie 
in the z — y plane. 

Figure [11] shows angle-slice projections of the gas sur- 
face density at times t — 23 Myrs (left panel) and t = 59 
Myrs (right panel). The projections are done along the x- 
axis, as before. In this projection both the disc axis of ro- 
tation (2- axis) and the outflow axis of symmetry are in the 
plane of the figure. It is apparent from Figure [11] that the 
outflow drives strong non-spherically symmetric motions in 
the infalling shell, sending gas on different non-circular or- 
bits inclined to the z-axis at various angles. 

The evolution of the shell is different from anything seen 
in the previous tests. Initially (see the right panel of the fig- 
ure) regions of the shell directly exposed to the feedback 
are compressed into a thin shell, which is slowly driven out- 
wards. However, because the shell is rotating, gas that was 
not previously exposed to feedback continuously replenishes 
the outflow's cone, and most of the shell eventually collapses 
inwards rather than being driven outwards, as is seen in the 
right panel of Figure 1111 

At the same time, the segments of the shell not exposed 
to feedback at all - regions along the diagonal line from the 
left upper corner to the right lower corner of the figure - 
approach the black hole initially almost freely falling (the 
right panel). At later times their radial motion is stopped 
by the centrifugal force. Some parts of the gas then appears 
to be pushed into the feedback cone by the ram pressure of 



the continuous inflow. These parts are blown away into the 
outflow. 

This geometrically complicated interplay of feedback 
and rotation results in a very strongly warped rotating struc- 
ture visible in Figure [12] at time t = 100 Myrs. Here we 
show two projections, one along the x-axis (left panel) and 
the other along the z-axis (right panel) for clarity. Note the 
circularising motions of gas at larger radii, settling into the 
disc. These regions of the disc are shielded from the black 
hole outflow by the inner disc regions. The inner disc regions 
are highly disturbed, however, with tangential, outflowing 
and inflowing motions present (see the central part of the 
right panel of Figure [12] ) . 

Figure [13] shows two later stages in the evolution of 
the system. The left panel of the figure is for time t — 160 
Myrs, whereas the right panel shows t = 400 Myrs. Appar- 
ently the feedback is eventually able to clear the directions 
in which it is operating to enforce a strong outflow in those 
directions. In short, parts of the disc perpendicular to the 
outflow do survive, whereas the less fortunate parts, in the 
path of the black hole wind, are blown away by the wind. 
This occurs on roughly the dynamical time of the original 
shell (Rout/cr « 190 Myrs). However this is probably a func- 
tion of the assumed disc temperature. A cooler and thicker 
disc could be more difficult to affect; more simulations of 
this are needed in the future. 

The resulting "truce" between the black hole outflow 
and the inertia of the disc is not an easy or natural one. 
The disc is warped and is still an evolving structure, with 
gas orbits in the inner part being nearly circular while those 
in the outer parts are more eccentric. The outflow also has 
a faster and slower parts, with the denser regions moving 
more slowly. Not all parts of the outflow reach escape ve- 
locities before becoming shielded, and some falls back along 
directions not exposed to the black hole outflow. 

Comparing this simulation to the one with isotropic 
black hole wind f 34.2[) , it is apparent that collimated feed- 
back leads the separation of the gas into disc and outflow re- 
gions to occur much later in time. Furthermore, this separa- 
tion is nowhere near as clear cut as in 34.21 In addition, while 
no gas was able to reach the accretion radius of R acc = 0.5 
kpc in the spherically symmetric feedback simulation, the 
misaligned feedback run resulted in about 1% of gas cross- 
ing the inner boundary. This value is still significantly less 
than the 4% of gas accreted in the control simulation with 
no feedback f 34.1[) . 

5.2 Outflow with a rotating axis of symmetry 

The final simulation that we present here has a set-up simi- 
lar to the previous simulation of misaligned collimated feed- 
back, except in this case the feedback axis is rotating around 
the x-axis according to: 

6Ut = Houtt , (14) 

where fl~^ t = 30 Myrs (the simulation in JET] corresponds to 
a constant ou t — tt/4). Thus the outflow axis becomes per- 
pendicular to its original position in 15 Myrs. Realistically, 
we might expect changes in the direction of black hole ro- 
tatio n to occur on even shorter timescales ([King fe Pringld 
l2006h , and so this simulation is probably still less complex 
than realistic bulges. 
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Figure 11. Angle-slice projected densities for the misaligned simulation at times t = 23 (left panel) and t = 59 (right panel) Myrs. 



Figure [14] shows the angle-slice projections of the sim- 
ulation at times = 47 (left panel) and t = 83 (right panel) 
Myrs. Figure [15] shows the same but at times = 120 (left 
panel) and t = 165 (right panel). There are certain simi- 
larities in gas dynamics between the present simulation and 
that with the misaligned collimated outflow (jjHJ), but there 
are also significant differences. One similarity is the pres- 
ence of gas inflow along directions not presently exposed to 
the black hole outflow, and outflow along some directions. 
However, there does not appear to be a settled structure in 
this case at all. Very strong velocity gradients and transient 
shocks occur in the present simulation. Such time dependent 
shocks may be efficient drivers of gas turbulence. 



6 FEEDBACK AND TURBULENCE 

Figure [16] shows the radial velocity structure for the SPH 
particles in the isotropic feedback simulation in which the 
shell had some angular momentum at a time of t = 236 
Myrs. As we remarked in M4.21 the feedback and the disc 
achieve a natural division of the available space on the 
spheres of influence. Along the directions close to the pole, 
feedback dominates. In the midplane, most of the disc par- 
ticles are shielded by the inner edge of the disc. The inter- 
mediate region is taken up by the wind driven off from the 
inner edge of the disc by the feedback. These features are 
clearly observable in Figure [16] in terms of the radial velocity 
component. 

We now consider the radial velocity plots for the simu- 
lation with a misaligned (stationary axis) outflow, fc|5.1l Fig- 
ure [17] shows that initially the radial velocity distribution is 
quite complex, but this settles into the warped disc - outflow 
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Figure 16. Radial velocity of SPH particles versus their radial 
distance to the black hole for the isotropic feedback simulation 
with rotation, see A4.2I and figure [10] The dominant feature is 
rotationally supported disc, and the points above it are the SPH 
particles in the outflow. 

structure fairly quickly already by the time t = 59 (lower 
panel) of the figure. With time this configuration evolves 
into a progressively more quasi-stationary state. 

In contrast, when the outflow is rotating ( Figure P"8]) . no 
clear cut steady state is reached. The radial velocity pattern 
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Figure 12. Left panel - As in Figure [TT] but for time t = 105 Myrs; note that the scale has changed. Right panel - same time but the 
density projection is along the z-axis. 
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Figure 13. As in Figure [TTl but for later times: t = 160 Myrs and t = 400 Myrs for the left and the right panels respectively. The outflow 
eventually evacuates the directions along which it acts; only the inclined part of the disc that is shielded from the feedback survives to 
late times. 
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Figure 14. Simulations with a rotating "jet" 


outflow. Left panel 


is for time t = 47 Myrs, and right is 


for t = 83 Myrs. 
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Figure 15. Same as Figure [14] but for later times: t = 120 Myrs and t = 165 Myrs for the left and the right panels, respectively. 
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Figure 17. Same as Figure fT6l but for simulation with an inclined 
collimated outflow f ^5.ip at two different times, t = 23 (upper 
panel) and t = 59 (lower panel) Myrs. See Figure [TT] for the edge- 
on view of the simulation. 



Figure 18. Same as Figure [TH but for simulation with a rotating 
collimated outflow ( ^5.2|) at two different times, t = 141 (upper 
panel) and t = 165 (lower panel) Myrs. The vertical feature in 
the lower panel is a strong transient shock. 



keeps evolving on time scales comparable to the rotation 
period of the jet at relatively late times. One observes de- 
velopment of occasional transient shocks such as the vertical 
feature in the bottom panel of Figure [18] The SMBH feed- 
back thus provides momentum input into complicated dif- 
ferential gas motions, probably pumping turbulence in the 
bulge. 

We believe this is in fact another way in which the 
SMBH can be affecting its host galaxy. The current consen- 
sus in the field of star formation points to the key role of tur- 
bulence in providing support against gravitational collapse 
of gas on large scales while allowing star formation to pro- 
cee d on smaller scales in dense filaments (for recent rev iews 
see iMac Low &; Klessenll2003 : iMcKee fe Ostriker]l200/t ). As 
turbulent motions are not equivalent to an isotropic pressure 



support (jDobbs et alJ l2QQ5h . the main effect of the turbu- 
lence is to delay or slow down star formation. While future 
higher resolution studies are needed to confirm the point, 
it appears intuitively clear that AGN-driven turbulence in 
the bulge can also hinder star formation in its host bulge. 
In this case the host is affected by the SMBH feedback in 
two ways: not only the gas is being blown away but also the 
fraction of gas turned into stars is reduced. 
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7 DISCUSSION 



7.1 Spherically symmetric models: successful 
feedback 

Spherically symmetric models can be considered a success 
to a certain degree since earlier analytical work is confirmed 
by the more detailed calculations here. First of all, t he an- 



alytical esti mate of the critical M CT mass obtained bv iKin 
(2003, 2005) agrees well with our numerical simulations, es- 
pecially for an initially static shell (see fc|2. 2|) . If the shell is 
already falling in, then the added inwardly directed inertia 
of the shell raises the M a mass by a factor of a few at most 
(g23] and g22}. 

The complications arising from a fin i te syste m size were 
recently pointed out bv lNavakshin et all (|2QQ9bh who noted 
that dynamical time in the bulge, Rb/cr, where Rb is the 
scale radius of the bulge, plays a key role in what happens. 
Black holes growing at the Eddington rate, MEdd, e-fold 
their mass on the Salpeter time, ts = MsH/MEdd ~ 4.5 x 
10 7 . Observationally, 



tdy 



= 17Myrs 



150 km s" 



(15) 



Thus "large" bulges, a ^ 200 km s _1 , have dynamical times 
comparable or larger than ts . Under the assumption that the 
duration of the black hole feeding event is comparable to the 
dynamic time of the bulge then, these "large" bulges allow 
sufficient time for their SMBH to grow significantly. These 
bulges would then have "fully grown" black holes saturated 
at their M a mass. On the other hand, lower mass bulges 
with velocity dispersions of a £ 100 km s _1 would evolve too 
rapidly for their SMBHs to "catch up". Therefore feedback 
from these SMBH would be overpowered by the inflow of 
gas from the bulge. The central regions of these lower mass 
bulges could be sites of copious star formation, leading to the 
birth of nucl ear star clusters. The latt er then reach their M a * 
masses (see iNavakshin et al.ll200 9b), replacing the SMBH 
black holes as the dominant objects in the bulges. 

Our simulations (*J3J) confirm this effect. The numerical 
experiment in ^3.11 has shown that black hole immersed in 
a 40 kpc bulge is able to grow enough to expel the gaseous 
shell when started from initial mass of M bh = 10 8 M (Fig 
[6]). At the same time, a black hole that is initially twice as 
massive in a 10 kpc bulge is unable to expel the shell ( fc|3.21 
Fig. [7J) that is less massive than that in £13.11 

One complication not considered by analytical models 
yet is that if the mass fraction of gas is significant, i.e. larger 
than the universal cosmological fraction of f g ~ 0.16, then 
the self-gravity of the gas can dominate the potential. The 
weight argument of course applies to self- gravitating shells 
as well, and they too can be unbound by a sufficiently strong 
black hole wind (e.g., see Figure [6}. However in this case our 
models are inconclusive as they do not include star forma- 
tion. When gas is blown away, the velocity dispersion of the 
self- gravitating shell should decrease. A detailed calculation 
including star formation is needed to establish whether or 
not the black hole in the remaining bulge follows the M — a 
relation. 



7.2 Spherically symmetric models: inconsistent 
SMBH feeding 

There are at least two major worries about the general use- 
fulness of spherically symmetric models. Large scale cos- 
mological simulations usually produce very complicated gas 
density an d velocity flows on to the centres of dark matter 
halos (e.g.. iKeres et al.ll2005h and so it is not clear whether 
or not such spherically symmetric flows ever take place. 

There is also a conceptual difficulty with the spherically 
symmetric models. In a self- consistent model, where black 
hole feedback is linked to the accretion of gas onto the hole, 
one cannot have accretion and outflow at the same time. Yet 
both are needed. We need accretion to grow and power the 
SMBH, but we also need outflow to eventually curtail the 
SMBH growth and expel the gas from the galaxy. 

Having accretion before, e.g., while the SMBH has not 
reached the M a mass, and outflow later, after it has done 
so, is not a comfortable option. Indeed, to change inflow to 
outflow for most of the gas in a bulge, we require momentum 
input of 



Mgas^e: 



MbulgeO" 



(16) 



where M gas is the mass of the gas which we assumed to 
be of order the bulge mass, Mbuige, and the escape velocity 
v esc ~ cr. The minimum amount of mass that the black hole 
needs to accrete to produce this much momentum outflow 
is 

P 



AM min = 



(17) 



For fiducial numbers a = 150 km/sec and e = 0.1, we have 

AM min = 0.005M bu i ge , (18) 

which i s comparable to the observed Mbh masses for most 
bulges (jHaring Rixll2Q04h . Now, if this gaseous mass is 
accumulated somewhere near the SMBH, e.g., in a disc that 
is immune to the feedback so that it can power the SMBH, 
then the question is why could this mass not be much higher 
or much lower? 

In other words, if SMBH feeding is local to a small 
scale region (for example, inside the SMBH influence radius, 
which is typically between a pc and a few tens of pc), where 
enough material is stored to feed SMBH while most gas is 
driven away, the casual link between feedback and accretion 
is broken, and there does not appear a reason to have an 
Mbh — cr relation. Thus spherically symmetric models, while 
providing excellent fits to the observed relations, do not nat- 
urally explain how the SMBH gets its fuel, or how it knows 
to stop growing once the many kpc scale shell starts to be 
blown away. 



7.3 Non-spherical models: clear mode for feeding 
but not feedback 

The second type of simulations that we considered here do 
not have spherical symmetry for one or more reasons. SMBH 
feeding in these models can be achieved (1) by features 
denser than the mean, i.e., discs or filaments, or (2) through 
regions not exposed to feedback at all. In the former case, if 
the SMBH-feeding gas has a column density (as seen from 
the black hole) much higher than the mean for the bulge, 
then it is only weakly affected by the momentum feedback. 
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In both of these cases there can be inflow to feed the SMBH 
and outflow to expel most of the gas in the bulge simul- 
taneously. However, if these dense features feed the SMBH 
despite it producing feedback at the maximum rate, then it 
is not clear why there should be a feedback mediated SMBH 
- galaxy link. 

This problem is already apparent in the simplest non- 
spherical case - the simulation of a rotating initially spher- 
ical shell of gas (S2J- In that simulation gas separates out 
in two well defined regions - the polar region where black 
hole outflow dominates and the disc midplane where gas is 
shielded from the black hole outflow by the inner edge of 
the disc. The edge of the disc is slowly "evaporating" into 
an outflow. The rest of the disc is hardly affected, and the 
orbits there are nearly circular. 

We have also set a numerical experiment in £15.11 where 
feedback outflow is confined to a broad cone misaligned with 
the rotation axis of the shell. The shell initial condition was 
identical to that explored in 34.21 We found, not surpris- 
ingly, that the resulting pattern of gas dynamics is even more 
complex than that of an isotropic feedback with a rotating 
shell. The inflow persists along directions in which black hole 
winds are not emitted, and the outflow prevails inside the 
feedback cone. The collimated rotating black hole outflow 
( N5.2|) is even more complex geometrically, with black hole 
feedback probably driving turbulence into bulge gas. 

7.4 Shortcomings and comparison to other 
numerical work 

There is a substantial body of work on numerical simula- 
tions of SMBH feed b ack in cosmological simulations, e.g., 
ISpringel et all (l2005h : lDi Matteo et all (|2005h : ISiiacki et al.1 
(2007). This paper complements this work. While address- 
ing similar issues, the scales involved in studying SMBH 
feedback in a cosmological context are very different. Cos- 
mological simulations have clear strengths - they naturally 
account for hierarchical build-up of galaxies assuming real- 
istic initial conditions and they include recipes for physical 
processes that we expect to be important for galaxy forma- 
tion such as cooling, star formation and feedback. There- 
fore they can provide a powerful tool with which to explore 
SMBH feedback. In comparison, our simulations are more 
modest in terms of the number of physical processes studied, 
making use of idealised initial conditions in static analytical 
potentials. 

On the other hand, the all-inclusive nature of cosmo- 
logical simulations mean that "individual" processes, such 
as black hole accretion and feedback, cannot be modelled in 
sufficient physical detail and require sub-grid physical mod- 
els to be assumed. As far as w e are aware, all the curr e nt cos- 
molog ical simulations rely on iBondi Hovlel (|l944h : iBondl 
formulation for accretion of gas. The latter formula- 
tion is physically inconsistent if accreting gas has angular 
momentum as gas simply cannot reach the SMBH due to 
the centrifugal barrier. Furthermore, at least the "quasar" 
mode of feedback is done in a very sketchy manner, with 
a fixed fraction of energy released by the quasar passed to 
the SMBH neighbour particles. The right fraction is found 
a posteriori to fit the observed Mbh — & relation and other 
observational constraints. 

Therefore, modest "small scale" simulations such as 



ours, especially when extended to include more physics and 
to reach in as close as possible to the SMBH may offer an- 
other important route to solving the mystery of SMBH and 
galaxy growth. Ideally, the small scale and cosmological sim- 
ulations should be merged to some degree in the future to 
allow as self-consistent a study of SMBH feedback as possi- 
ble. 



8 CONCLUSIONS 

Our main conclusions are 

• Predictions of spherically symmetric models with black 
hole feed back t i ed to Eddington limit luminosity as in the 
models of King (2003, 2005) are confirmed numerically. Self- 
gravity of the gas complicates evolution of the system, and a 
self- consistent treatment of star formation and its feedback 
is necessary for further progress. 

• As suggested earlier based on analytical arguments, it 
is the dynamical time in the bulge, Rb/cr, that determines 
whether or not the SMBH can reach its limiting M a value. 
Central regions of smaller bulges, where the dynamical time 
is observed to be shorter than the Salpeter time, could be 
smothered with infalling gas despite the SMBH feedback. 
This process may be the origin of the nuclear star clusters 
in "smaller" galaxies. 

• A net angular momentum in the shell is essential in 
determining the fate of the shell and the SMBH feeding. 
There is no well defined M a mass in that case since the mo- 
mentum thrust required is different in different directions. 
Gas near the symmetry axis is blown out easier than in the 
spherically symmetric case, whereas gas settled into a disc 
requires ~ R/H more thrust than in the latter case. 

• If the SMBH is fed through a several kpc scale disc, the 
SMBH mass would have to be ~ R/H larger than in the 
spherical case to expel the feedback-resistant self- shiel ding 
disc. T his directly contradicts the recent observations of iHul 
(|2009h that show that pseudo-bulges have lighter SMBHs 
than their classical counter-parts. This may imply that black 
holes are not fed by large (kpc or larger) discs but rather by 
flows with a small specific angular momentum. 

• We also noted that SMBH outflows in a realistic, i.e., 
aspherical, situation may pump turbulence (differential mo- 
tions) in the bulge. Turbulence is known to hinder star for- 
mation. Therefore we believe that AGN feedback not only 
expels the gas from the galaxy but it may also reduce the 
amount of mass turned into stars during bulge formation. 
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